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\D '. Abstract 

The K^K^ scattering length is calculated in fully-dynamical lattice QCD with domain-wall va- 
^ _ lence quarks on the MILC asqtad-improved gauge configurations with rooted staggered sea quarks. 

' Three-flavor mixed-action chiral perturbation theory at next-to-leading order, which includes the 

■ leading eff'ects of the finite lattice spacing, is used to extrapolate the results of the lattice calculation 

to the physical value of m^+ / fj^+. We find m^+ ax+K+ = —0.352 it 0.016, where the statistical 



X 



and systematic errors have been combined in quadrature. 
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I. INTRODUCTION 



Strange hadrons may play a crucial role in the properties and evolution of nuclear material 
under extreme conditions [l|. The interior of neutron stars provide one such environment in 
which the densities are high enough that it may be energetically favorable to have strange 
baryons present in significant quantities, depending upon their interactions with non-strange 
hadrons. Further, it may be the case that a kaon condensate forms due to strong interactions 
between kaons and nucleons Unfortunately, the theoretical analysis of both scenarios 
is somewhat plagued by the limited knowledge of the interactions of strange hadrons with 
themselves and with non-strange hadrons. 

Heavy- ion collisions, such as those at the BNL Relativistic Heavy Ion Collider (RHIC), 
also produce nuclear material in an extreme condition. Recent observations suggesting the 
formation of a low- viscosity fluid are quite exciting as they provide a first glimpse of matter 
not seen previously. The late-time evolution of such a collision requires an understanding 
of the interaction between many species of hadrons, not just those of the initial state, 
including the interactions between strange mesons and baryons. While pion interferometry 
in heavy-ion collisions is a well-established tool for studying the collision region (for recent 
theoretical progress, see Refs. [1, 0, lEI), the STAR collaboration has recently published 
the first observation of neutral kaon {K'^) interferometry [6|. In the analysis of K^-K^ 
interferometry, the non-resonant contributions to the final state interactions between the 
kaons were estimated using three-flavor {SU{?>)l ® SU{3)ji) chiral perturbation theory 
PT), the low-energy effective field theory of QCD. Given the sometimes poor convergence 
of SU{3)l ® SU{3)ji x-PT due to the relatively large kaon mass compared to the scale of 
chiral symmetry breaking (A^ ~ 1 GeV), particularly in the baryon sector, it is important 
to be able to verify that the non-resonant contributions to f^i^-scattering are indeed small, 
as estimated in X"PT. 

In this work we present the first lattice QCD calculation of the K~^K^ scattering length. 
The calculations are performed on the coarse MILC lattices (with a preliminary calculation 
on one ensemble of the fine MILC lattices) and three-flavor mixed-action x-PT (MA^- 
PT), which includes the leading-order lattice-spacing effects, is used to extrapolate to the 
continuum and to the (isospin-symmetric) physical value of the meson masses. We find that 
at the physical value of mx+ / /x+ 

mK+ aK+K+ = -0.352 ± 0.016 , (1) 

where the statistical and systematic errors have been combined in quadrature. 

The TT, K and t] are identified as the pseudo-Goldstone bosons associated with the spon- 
taneous breaking of the approximate chiral symmetry of quantum chromodynamics (QCD), 
and therefore the form of their interactions is highly constrained. In fact, at leading order 
(LO) in X"PT, the scattering of two of these mesons is uniquely determined 0]. Corrections 
to the LO scattering amplitude arise in a systematic expansion about the chiral limit [sjo], 
scaling generically as (m^^^/A^)*^ where n counts the order in the chiral expansion 10|. 
For obvious reasons, SU (3) x-PT is expected to converge more slowly than two-flavor x-PT. 

There is a wealth of phenomenological and theoretical knowledge concerning low-energy 
TTTT scattering. The chiral extrapolation formulae for tttt scattering are known to two loops, 
or next-to-next-to-leading order (NNLO_K in both SU{2) [l| and SU{3) [Ts) x-PT. 



Combined with a Roy equation analysis 114,1 llSl . Il6| , this has allowed for remarkably-precise 



determinations of the two s-wave vrvr scattering lengths [17|, Il8|, Il9| . In the case of the Kn 
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systems, the extrapolation formulae for the scattering amplitudes are known to one [201 . 
21 . 22 | and two loops [23| and have allowed for theoretical predictions of the / = 1/2 
and 7 = 3/2 scattering lengths. However, the uncertainty in these theoretical predictions 
is substantial. There are proposed experiments to study the Kn atoms by the DIRAC 
collaboration [13] at CERN, J-PARC and GSI, which will significantly reduce the uncertainty 
in these scattering lengths. To date, there have been no experimental determinations of the 
7 = 1 7^7^ scattering length, aj^, but recently it has been calculated at next-to- leading 
order (NLO) in x-PT 

The methods for studying two-particle interactions in a finite Euclidean volume are well 
known (26l . 27 . [28l . 29| . The interaction energy of two hadrons in a finite volume uniquely 
determines pcot5(p), and hence their scattering amplitude, below kinematic thresholds. 
The scattering parameters, such as the scattering length and effective range, can then be 
determined from calculations of pcot5(p) over a range of energies. These methods paved 



the way for 
decade ago | 



)ioneering quenched QCD calculations of two-particle interactions a little over a 
|33|, |3J]. Since then, there have been a number of additional quenched 

Si. The first 



36, 37 



39, m, 41 



calculations of the 7 = 2 tttt scattering length [3| 
dynamical calculation of vrvr interactions (including the phase-shift) was performe d by the 
CP-PACS collaboration with two flavors (n/ = 2) of improved Wilson fermions [43| and 
pion masses in the range 500 < m^r < 1100 MeV. Recently, dynamical calculations of the 
7 = 2 TTTT scattering length with three flavors of light quarks (n/ = 2-1-1) were performed 
with pion masses in the range 300 < < 500 MeV. These calculations used a mixed- 
action scheme of domain-wall valence fermions on asqtad-improved staggered sea fermions 



at a single lattice spacing of 6 ~ 0.125 fm [44, |45(], and used mixed- action X"PT (MAy-PT) 
(which describes the finite-lattice spacing effects) calculations of the scattering length [46| 
to extrapolate to the physical meson masses. Only recently has the first calculation of the 
7 = 3/2 7^77 scattering length in quenched QCD been performed and a fully-dynamical 
Hf = 2 + 1 calculation (48| followed shortly afterward. When combined with X"PT, this 
latter calculation allowed for a simultaneous prediction of the 7 = 1/2 and 7 = 3/2 Kir 
scattering lengths using the NLO extrapolation formulae j48| . 

This paper is organized as follows. Section [ITl contains the details of our mixed-action 
lattice QCD calculation. Discussion of the relevant correlation functions and an outline 
of the methodology and fitting procedures can also be found in this section. The results 
of the lattice calculation and the analysis with MA^-PT are presented in Section IIIII In 
this section, the various sources of systematic uncertainty are identified and quantified. In 
Section [IV] we conclude. 



II. METHODOLOGY AND DETAILS OF THE LATTICE CALCULATION 



In calculating the K^K^ scattering length, the mixed-action lattice QCD scheme developed 
by the LHP Collaboration (49l . 50 was used in which domain- wall quark jsil . 52, 53, 
propagators are generated from a smeared source on = 2 -t- 1 asqtad-improved [56|, i57| 
rooted staggered sea quarks j58|. To improve the chiral symmetry properties of the domain- 



60, iMi was used in the gauge links of 



wall quarks, hypercubic-smearing (HYP-smearing) [59l . 
the valence-quark action. In the sea-quark sector, there has been significant debate regarding 
the validity of taking the fourth root of the staggered fermion determinant at finite lattice 



spacing [62, l63|, l6J, l65|, l66|, l67|, l68|, l69|, [TOj, LDJ, l72|, [73|, [7^, l75| . While there is no proof, there 



are arguments to suggest that taking the fourth root of the fermion determinant recovers 
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TABLE I: The parameters of the MILC gauge configurations and domain-wall propagators used in 
this work. The subscript / denotes light quark (up and down), and s denotes the strange quark. The 
superscript dwf denotes the bare-quark mass for the domain-wall fermion propagator calculation. 
The last column is the number of configurations times the number of sources per configuration. 



Ensemble 


bmi 


bnis 


T dw f 

brrij^ ' 


T dw f 

bms 


10'^ X bnires 


# of propagators 


2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


1.604 


468 


X 


16 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


1.552 


658 


X 


20 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


1.239 


486 


X 


24 


2064f21b681m030m050 


0.030 


0.050 


0.0478 


0.081 


0.982 


564 


X 


8 


2896f2b709m0062m031 


0.0062 


0.031 


0.0080 


0.0423 


~ 0.25 ^ 


506 


X 


1 



"Computed by the LHP collaboration. 
'Estimated on a small number of configurations. 



the contribution from a single Dirac fermion ^. The results of this paper assume that the 
fourth- root trick recovers the correct continuum limit of QCD. 

The present calculations were performed predominantly with the coarse MILC lattices 
with a lattice spacing of 6 ~ 0.125 fm, and a spatial extent of L ~ 2.5 fm. On these 
configurations, the strange quark was held fixed near its physical value while the deg enerate 



light quarks were varied over a range of masses; see Table [T] and Refs. |77|, l78|, l79|, l80|, ISl 
for details. Further, preliminary calculations were performed on 506 configurations of one 
fine MILC ensemble. On the coarse MILC lattices, Dirichlet boundary conditions were 
implemented to reduce the original time extent of 64 down to 32 and thus save a factor of 
two in computational time. While this procedure leads to minimal degradation of a nucleon 
signal, it does limit the number of time slices available for fitting meson properties. By 
contrast, on the fine MILC ensemble, anti-periodic boundary conditions were implemented 
and all time slices are available. To determine the light-quark masses, the domain-wall pion 



was tuned to the lightest staggered pion to within a few percent |49|, |50]. This choice is 
somewhat arbitrary as the partially-quenched and mixed-action effective field theories exist 
to describe this and other choices |82| (provided the meson masses remain in the chiral 
regime), with the expression for the 1 = 1 KK scattering length determined to NLO in 
X-PT, PQx-PT and MA^-PT in Ref. j25|.^ The choice of tuning to the lightest taste 
of staggered meson mass, as opposed to one of the other tastes, provides for the "most 
chiral" domain-wall mesons and therefore reduces the error in extrapolating to the physical 
point. The mass splitting between the do main- wall mesons and the staggered taste-identity 
mesons, which characterizes the unitarity violations present in the calculation, is then given 



by |8i,l85 



m^^-m^^^^^62Ai = 0.0769(22) (l.u.) 



0.0295(27) (l.u.) fine. (2) 



^ For a nice introduction to staggered fermions and the fourth-root trick, see Ref. 76 1. 

^ The PQx-PT and MAx-PT expressions for the 7=1 KK scattering length are identical in form at NLO. 
This is not unique to this quantity and can be understood on more general grounds, as mixed-action 
theories with chirally-symmctric valence fermions exhibit many universal features 83[. 
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In order to determine the interaction energy between the two kaons, both the single- 
kaon, CK+(t), and two-kaon, CK+K+{p,t), correlation functions were computed, where t is 
the Euchdean time separation between the hadronic source and sink operators and p denotes 
the magnitude of the equal and opposite spatial momentum of each kaon. The single-kaon 
correlation function is 

CMt) = J2{K'{t,^) K^{0,0)) , (3) 



X 



where the sum over all spatial sites projects onto the zero-momentum state, p = 0. A 
correlation function which projects onto the K^K^ s-wave state in the continuum limit is 

Ck^k^M = J2 $^e^'"^^-^)(i^-(t,x) K-{t,y) K+{0,0) A^+(0,0)). (4) 

|p|=P x,y 

In eqs. ([3]) and (jl]), K^{t,:x.) = x)75u(t, x) is a Gaussian-smeared interpolating field 
for the charged kaon. In the relatively-large spatial volumes used in the calculation, the 
interaction energy between the two kaons is a small fraction of the total energy, which is 
dominated by the kaon masses. To determine this energy, the ratio of correlation functions, 

was constructed, with the arrow denoting the large-time, infinite-number-of-gauge- 
configurations limit (far from the boundary). Due to the periodic boundary-conditions 
imposed on the propagators computed on the fine lattices, the correlation function 
become a single cosh function far from the source, while the K^K^ correlation function 
become the sum of two cosh's, one depending upon mK+ and the other depending upon 
Ek+k+, leading to a non-trivial form for Gx+K+iPii)- As an alternative method to calcu- 
lating the interaction energy (and a check of the systematics) , a Jackknife analysis of the 
difference between the energies extracted from the long-time behavior of the double- and 
single-kaon correlation functions individually was performed, finding results in agreement 
with those determined from eq. (jS]). The interaction energy is related to the two-particle 
energy eigenvalues and twice the kaon mass. 



AE„ = - 2mK = 2^pl + ml - 2mK • (6) 

In the absence of interactions, the energy levels occur at values of the momenta p = 27rj/L 
where j is an integer-triplet, corresponding to the allowed single-particle momentum modes 
in a cubic volume. In the interacting theory, the two-particle eigen- moment a, are shifted 
from these values and can be determined from eq. (jH]) and the calculated interaction energy. 
The Liischer formula [2^, 27, 28, 29 1 can then be used to determine the infinite- volume 
scattering parameters from the real part of the inverse scattering amplitude by solving the 
equation 

pcot.(rt^-ls(§), (7) 

which is valid below the inelastic threshold. The regulated three-dimensional sum is 86 

|j|<A 



S(^)-E^p3^-4-A, (8) 
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TABLE II: Masses, energies and scattering lengths determined from the lattice calculation. The 
first uncertainty assigned to each quantity is statistical, determined with the Jackknife procedure, 
and the second uncertainty is an estimated fitting systematic. 



Quantity 


mi = 0.007 mi = 0.010 mi = 0.020 mi = 0.030 


mi = 0.0062 


b tUt^ 
Fit Range 


0.1846(4)(2) 0.2226(4)(3) 0.3104(3)(15) 0.3747(4)(8) 
8-14 9-13 9-15 6-13 


0.1453(5)(13) 
17-39 


b rriK 
Fit Range 


0.3680(4)(4) 0.3776(3)(4) 0.4046(3)(13) 0.4300(4)(3) 
7-11 9-15 9-15 9-13 


0.2458(5)(13) 
20-34 


m-n/fK 

rriK/fK 


1.712(4)(3) 2.069(3)(5) 2.835(3)(11) 3.335(4)(9) 
3.412(5)(4) 3.509(3)(6) 3.695(3)(10) 3.827(4)(9) 


1.978(15)(12) 
3.344(19)(21) 


Fit Range 


0.00619(30)(32) 0.00663(15)(35) 0.00606(14)(22) 0.00613(19)(10) 
12-17 10-16 11-17 12-17 


0.00437(36)(105) 
18-34 


(MO) 


-0.448(19)(20) -0.497(10)(22) -0.523(10)(23) -0.590(15)(21) 


-0.391(28)(82) 


which runs over all triplets of integers j such that j < A and the limit A — oo is implicit. 
The scattering parameters are then related to p cot 6{p) through the effective-range expansion 

p cot 6{p) = ^ + ^r/ + C(/) , (9) 



where a is the scattering length and r is the effective range. For naturally-sized scattering 
lengths and small interaction momenta, p cot 6{p) is predominantly given by the inverse 
scattering length. 



III. ANALYSIS AND THE CHIRAL AND CONTINUUM EXTRAPOLATIONS 

It is convenient to present the results of our calculation with "effective scattering length" 
plots, determined from the ratio of correlation functions, 

and similarly on the fine ensemble. For each time slice, AE{t) is inserted into eq. ([7]) 
which yields a scattering length at each time slice, ax+K+it)- To remove any scale-setting 
ambiguities, the scattering length is multiplied by the "effective" kaon mass, mx{t). The 
effective scattering length plots associated with each lattice ensemble are shown in Fig. [21 
The statistical errors are determined from a Jackknife analysis, while the quoted systematic 
errors are estimated from both the range of fits as well as the two methods of determining 
the interaction energy described in Sec. [Tll In Table [III the calculated values of the meson 
masses, decay constants, two-particle energy shifts and scattering lengths are presented. 
Effective kaon mass plots and effective scattering length plots are shown in Fig. [1] and 
Fig. [21 respectively. 



6 



0.374 
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-^^ 0.370 
S 0.368 
^ 0.366 
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[i: 0.38 
W 0.378 

s 
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_Q 0.404 

0.402 

0.434 
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FIG. 1: The effective mx+{t) plots. The solid black lines and shaded regions are fits with 1-a 
statistical uncertainties (Table\l^. The dashed lines correspond to the statistical and systematic 
(Table\l^ uncertainties added in quadrature. 
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FIG. 2: The effective K^K^ scattering length times the effective as a function of time 

slice. The solid black lines and shaded regions are fits with 1-a statistical uncertainties (Table IT/j) . 
The dashed lines correspond to the statistical and systematic (Table uncertainties added in 
quadrature. 
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A. Mixed- Action x-PT at One Loop 



The lattice QCD calculations performed in this work are isospin symmetric, m„ = m^, and 
do not include electromagnetism. Therefore isospin is a good quantum number. Having 
computed the K^K^ scattering length at a number of unphysical pion masses and at a 
finite lattice-spacing, isospin-symmetric MA^-PT is used to extrapolate to the physical 
(isospin-symmetric) meson masses and to the continuum. 



In Ref. j25j, the expression for the 1 = 1 KK scattering length was determined to NLO 
in X"PT, including corrections due to mixed-action lattice artifacts. As with the 1 = 2 vrvr 



scattering length [46(] , it was demonstrated that when the mixed-action extrapolation formula 



is expressed in terms of the lattice-physical parameters computed on the lattice, ^ m^r, rriK 
and fx, there are no lattice-spacing-dependent counterterms at 0(&^m|-) or 

There are finite lattice-spacing-dependent corrections, proportional to b'^Ai, and therefore 
entirely determined to this order in MAx-PT. Again, as with the / = 2 vrvr system, the 
NLO MA formula for m^^aj^ does not depend upon the mixed valence-sea meson masses, 
and therefore does not require knowledge of the mixed-meson masses jSTj. This allows for a 
precise determination of the predicted MA corrections to the scattering length. At NLO in 
MAx-PT, the scattering length takes the form 



mxaKKib ^ 0) 



m 



K 



+ Cvln 



m 



K 



(4vr/ 



~2 



m 



X 



m. 



+ Css In 



m 



K 



+ Co- 32(47r)^L^t(/^) 



fill 



where the various coefficients, Ci, along with and m^^, can be found in Appendix E 
of Ref. To account for the predicted MA corrections, one can either use eq. ( TTT]) to 

directly fit the results of the lattice calculation (Table HT]) or one can determine the quantity 



MA 



rriKaKK 



XPT 



(12) 



collected in Table IIIII and Table llVt subtract this from the results of the lattice calculation 
and use the NLO X"PT expression for the scattering length. 



m 



K 



m 



K 



(4vr/. 



2 In 



m 



K 



+ 



2(20m2^ - llm^) 
27{m'i — ml) 



2ml 



\m. 



mt 



In 



ln(^j-y-32(4vr)^Lri(^) 



• (13) 



As there is only one counterterm at NLO, it can be determined on each ensemble. In order 
to carry out this analysis, further sources of systematic errors are identified; higher-order 
effects in the chiral expansion, [S.tsinlo{"^k(^^kk)'^ exponentially-suppressed finite-volume ef 
fects 
A 



11 CApciiioiuii, i-^]<j^i^Qynii^C^ ^ 

lS.py{mKO^KK)'^ residual chiral symmetry breaking effects from the domain- wall action. 



[mKa^KK) 



and the error in truncating the effective-range expansion with the inverse 



^ Quantities calculated directly from the correlation functions are denoted as lattice-physical parameters. 
These are not extrapolated to the continuum, to infinite volume nor to the physical quark mass point. 
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TABLE III: The continuum limit of the scattering length at the physical point on the coarse MILC 
lattices, the extracted counterterm that enters at NLO in x-PT, and the various systematic un- 
certainties that have been identified beyond those associated with fitting. The correction factors, 
Ai, are defined in the text. The first uncertainty associated with each scattering length is statis- 
tical, the second is the systematic uncertainty from Table O (ind the third is from the systematic 
uncertainties presented in this table (combined in quadrature). The first uncertainty associated 
with each L^xkil^ ~ f^) is statistical, while the second is systematic (all systematics combined in 
quadrature). 



Quantity 


mi = 0.007 


mi = 0.010 


mi = 0.020 


mi = 0.030 


Afv (iT^Ka^KK) 


-0.0067(14) 
±0.016 
±0.001 
±0.007 
± 0.008 


-0.0062(16) 
±0.019 
±0.001 
±0.006 
±0.008 


-0.0052(19) 
±0.028 
±0.000 
±0.005 
±0.008 


-0.0048(21) 
±0.037 
±0.000 
±0.004 
±0.007 


mK+aK+K+ 
(6-0) 


-0.441(19)(20)(19) 


-0.491(10)(22)(22) 


-0.518(10)(23)(30) 


-0.585(15)(21)(38) 


32(4^)2l^=^(/^) 


7.3(5)(8) 


6.8(3)(8) 


7.7(2)(8) 


7.4(3)(8) 



TABLE IV: The continuum limit of the scattering length at the physical point on the fine MILC 
lattices, the extracted counterterm that enters at NLO in x-PT, and the various systematic uncer- 
tainties that have been identified beyond those associated with fitting. The correction factors and 
uncertainties are discussed in the caption of Table [77/1 



Quantity 


mi = 0.0062 


Ama {mKa^KK:) 


-0.0048(15) 




±0.013 


A-FV ("^xajfi) 


±0.001 


Amres {mKa^^K) 


±0.004 


Arange {^K^KK^ 


±0.004 


mK+aK+K+ 


-0.387(28)(82)(14) 


(6-0) 




32(47r)2L^=^(/^) 


8.4(9)(2.6) 



scattering length, ArangeimRdKK)- These various sources of systematic uncertainty, as well 
as the predicted mixed- action corrections, the adjusted scattering lengths and the deter- 
mined values of L^]^(/i) are given in Table UTTl and Table HVl In the following sections, each 
source of systematic uncertainty is addressed in turn. 

1. NNLO x-PT Corrections 

The NNLO extrapolation formula for rnxdj^x does not exist, and therefore estimates of 
contributions from higher order in the chiral expansion are limited to power-counting argu- 
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ments. A conservative estimate is provided by 

^NNLO {^kO-^Kk) = 



(47r/. 



In 



m 



K 



K 



(14) 



and the resulting uncertainties are given in Table IIIII and Table IIVI 



2. Finite-Volume Effects in Mixed- Action x-PT 

Liischer's relation between the two-particle energy levels in a finite volume and their infinite- 
volume scattering parameters receive exponential corrections which depend upon the lattice 
size and the lightest particle in the spectrum, and generically scale as e~™'"^. In Ref. 88| . 
the exponential volume corrections were determined for the J = 2 tttt system. Using these 
methods, one can also determine the exponential volume corrections for the mixed-action 
K^K^ system. These exponentially-suppressed volume corrections are formally sub-leading 
compared to the effective-range corrections which have not been included, and provide an 
estimate of the finite-volume corrections. These terms are denoted as 



FV 



rriKa^KK 



ooV 



(15) 



and are collected in Table IIIII and Table IIVI 



3. Residual Chiral Symmetry Breaking 

The NLO mixed- action formula, eq. (|TT|) . as well as the corrections of Table UTTI and Table HVl 
were derived assuming valence fermions with perfect chiral symmetry. However, domain-wall 
fermions are necessarily implemented with a finite fifth-dimension which induces residual 
chiral symmetry breaking. The leading contributions from this residual chiral symmetry 



breaking can be parameterized with a residual quark mass |53|, |54 



dw f dw f I 

rrii ^ rrii + rrires ■, 



m 



mf^ + mres . (16) 



However, by expressing the MA^-PT formula in terms of the lattice-physical meson masses, 
the dominant contribution from these m^es terms are automatically included. This leaves 
corrections at NLO (assuming m^es ~ mq in the expansion), some of which have undeter- 
mined coefficients. Naive dimensional analysis [89| can be used to estimate the size of these 
terms, 

A,„..(m,a,,)=±^^-^--, (17) 
which are shown in Table IIIII and Table IIVI 



^. Range Corrections 

When the spatial dimensions of the lattice are large compared to the range of the interaction, 
and the scattering length is of natural size, as is the case for K^K^ scattering at the quark 
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FIG. 3: mx+CLx+K+ versus fnj^+ / fj^+ . The points with error-bars are the results of this lattice 
calculation (not extrapolated to the continuum) on both the coarse and fine MILC lattices. The 
solid curve corresponds to the tree-level prediction ofx-PT, and the point denoted by a star and its 
associated uncertainty is the value extrapolated to the physical meson masses and to the continuum. 
The smaller uncertainty associated with each point is statistical, while the larger uncertainty is the 
statistical and fitting systematic combined in quadrature. 



masses used in this work, the effective range first enters at 0{L~^) in the expansion of the 
two-hadron energy in powers of 1/L. Therefore, neglecting the effective-range parameter 
introduces a ~ 0.2% uncertainty in the extracted values of m^+a^+^+, assuming r ~ 
l/(2m^) for nin <C mx- To be conservative, a 1% systematic uncertainty due to the neglect 
of the effective range is assigned to the scattering length determined on each ensemble. 



B. Extrapolation to the Physical Point 

Calculations on the four coarse lattice ensembles yield pion and kaon masses of approximately 
(m^,mi^) ~ (290,580), (350,595), (490,640) and (590,675) MeV. The chiral expansion will 
converge better for smaller meson masses, and one method to examine the convergence of 



the chiral expansion is to selectively "prune" the heaviest data sets |44l . I45l . I48l |. This is 
done by first determining Lj^(/i = fx) by fitting to all four data points (fit A), then 
removing the heaviest point and fitting (fit B) and finally removing the heaviest two points 
and fitting (fit C). The results of these fits are collected in Table IVl The extracted values 
of 

^KK from each of the fits are consistent with each other within the uncertainties. In 
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TABLE V: The results of fitting three-flavor MAx-PT at NLO to the computed scattering lengths, 
as described in the text. The values of mj^+ax+K+ o,'^^ those extrapolated to the physical (isospin- 
symmetric) meson masses and to the continuum. The first uncertainty is statistical and the second 
is systematic (as described in the text). 



FIT 


32(4vr)2Lri(/i.) 


mK+aK+K+ (extrapolated) 


xVdof 


A 


7.3(1)(4) 


-0.347 ± 0.003 ± 0.009 


0.22 


B 


7.3(2)(5) 


-0.347 ±0.004 ±0.011 


0.32 


C 


6.9(2)(6) 


-0.355 ± 0.005 ± 0.013 


0.14 



analogy with the comparison convention employed for 7r+7r+, the lattice data is extrapolated 
to the physical values of = 0.8731 ± 0.0096, rriK+l fK+ = 3.088 ± 0.018 and 

^■nl fK+ = 3.425 ±0.0019 assuming isospin symmetry, and the absence of electromagnetism. 
Taking the range of values of L^j^^ spanned by these fits, we find 

mK+aK+K+ = -0.352 ±0.016 , 32{ATrfLj^^{fj. = fx) = 7.1 ±0.7, (18) 

where the statistical and systematic errors have been combined in quadrature. The results 
are shown in Fig. [31 It is somewhat surprising that the calculated scattering lengths are 
consistent, within uncertainties, with tree-level X"PT. This was also found to be the case for 
TT+TT"'' scattering even at large pion masses. 



C. Comparing K^K^ scattering with 7r+7r+ scattering 



A comparison between the lattice calculations of tt^tt"'" [45j and K^K^ scattering lengths 
allows for a study of fiavor-Sf/ (3) breaking in the scattering amplitude due to terms that are 
beyond NLO in X"PT. The linear combination of Gasser-Leutwyler coefficients contributing 
to the / = 2 TTTT scattering length at NLO is the same as the combination contributing to 



the 1 = 1 KK scattering length [25 



L'^=^{fi) = Lifi^^) . (19) 

To compare extractions of these counterterms, the scales at which they are evaluated must 
be the same, and the scale-dependence of Lj^(/i) is 

32(47r)2Lri(/.) = 32(47r)^Lri(/.o) - f In (^Q . (20) 

The counterterms extracted from the mixed-action lattice calculations are shown in fig. H] 
as a function of m^//^. It is clear that while there appears to be a difference between 
^K^xifn) ^-nd L^^^(/^), more precise calculations of both scattering lengths, particularly at 
the lightest pion masses, are required for further exploration of higher order terms in the 
chiral expansion. 



IV. DISCUSSION 

We have presented results of a lattice QCD calculation of the K^K^ scattering length 
performed with domain-wall valence quarks on asqtad-improved MILC configurations with 
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10 



9 




m /f 



FIG. 4: 32(47r)^L;^(/7r) (circles) and 32(47r)^L^^^(/^) (squares) versus m-,^/ f-,^. The sets have 
been slightly displaced horizontally for convenience of viewing. The smaller uncertainty associated 
with each point is statistical, while the larger uncertainty is the statistical and fitting systematic 
combined in quadrature. 



2+1 dynamical staggered quarks. The calculations were performed on the coarse MILC 
lattices with a lattice spacing of 6 ~ 0.125 fm (with a preliminary calculation on one ensemble 
of the fine MILC lattices with h ~ 0.09 fm) and at a single lattice spatial size of L ~ 2.5 fm. 
One-loop MAx-PT with three flavors of light quarks was used to perform the chiral and 
continuum extrapolations. Our prediction for the physical value of the K^K^ scattering 
length is mK+ clk+k+ = —0.352 ± 0.016, and we emphasize once again that this result 
rests on the assumption that the fourth-root trick recovers the correct continuum limit of 
QCD. Deviations from Weinberg's tree-level prediction are found to be surprisingly small, 
consistent with the lattice calculations of the 7t~^7:~^ scattering length at heavier pion masses. 
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